The fate of ID spin-charge separation away from Fermi points 
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We consider the dynamic response functions of interacting one dimensional spin- 1/2 fermions at 
arbitrary momenta. We build a nonperturbative zero-temperature theory of the threshold singulari- 
ties using mobile impurity Hamiltonians. The interaction induced low-energy spin-charge separation 
and power-law threshold singularities survive away from Fermi points. We express the threshold 
exponents in terms of the spinon spectrum. 
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The low-energy excitations of interacting spin- 1/2 
fermions confined to one dimension (ID) is well repre- 
sented by two collective bosonic modes. These modes 
are the quantized waves of spin and charge densities. 
Their spectra are linear; the corresponding velocities, v s 
and Vr, differ from each other. A microscopic consider- 
ation [l], Q of the repulsive interaction between spinful 
fermions leads to v c > v s . 

The spectra of the collective modes can be probed in a 
momentum-resolved tunneling 0, Q] or in an ARPES [B| 
or photoemission [fjj experiment. In these methods, a 
spin-1/2 fermion with a given momentum tunnels into or 
out of the studied system. The tunneling inevitably per- 
turbs each of the two collective modes. For example, in 
the case of a low-energy particle (hp 3> k — kp > 0), the 
small difference k — kp is shared between the excitations 
of the two modes. The Luttinger liquid (LL) theory pre- 
dicts 0] that at given k the tunneling probability is sin- 
gular at energies v s (k—kp) and v c (k—kp), corresponding 
to the entire momentum k — kp given to the "spinon" or 
"holon" belonging to the spin and charge mode, respec- 
tively. The exponents of the two power-law singularities 
depend on a single number, the LL parameter K c for the 
charge mode. The two sharp peaks in the momentum- 
resolved tunneling probability at energies associated with 
excitation of the two modes, to — v SiC ■ (k — kp), are the 
hallmark of the spin-charge separation in the LL. 

The momentum-resolved tunneling rate is propor- 
tional to the fermionic spectral function A{k,uj) — 
iRe / dtdxe lu,t - lkx 0{t) ({^[x,t),^ [0,0)}). Recently, 
considerable progress was achieved in the analytical the- 
ory of dynamic responses of a ID system away from 
the Fermi points for spinless fermions @-[llj]. The de- 
veloped methods map the ID dynamic response prob- 
lem near the edge of support onto the "mobile quan- 



tum impurity" effective Hamiltonian 12|. For spinless 
fermions in the weak-interaction limit, one may consider 
the generic spectrum of free fermions exactly, while treat- 
ing their interaction perturbatively 0, Q- For example, 
at |fc| < kp the threshold coincides with the spectrum 
of a hole, which can be thought of as a mobile quantum 
impurity. Because of the interactions, it "shakes up" the 



fermions in the vicinity of Fermi points, leading to the 
orthogonality catastrophe and to the power-law behavior 
of A{k,w) at the threshold The perturbation the- 

ory allows one to identify the quantum numbers of the 
impurity and to match the phenomenological theory of 
threshold exponents flfj| valid at any interaction strength 
with the weak- interaction limit. 

In this Letter we build a nonperturbative zero- 
temperature theory of threshold singularities of dynamic 
responses of ID spin-1/2 fermions at arbitrary k, shed- 
ding light on the fate of the spin-charge separation away 
from Fermi points. The obvious difficulty of the prob- 
lem lies in the appearance of distinctly different spin and 
charge modes at any interaction strength. For weak in- 
teractions, these modes are degenerate, which renders 
perturbation theory inapplicable. We find that interac- 
tion induced spin-charge separation survives away from 
the Fermi points, and manifests itself in the choice of 
quantum numbers of the mobile impurity: for repulsive 
interactions, it only carries spin but no charge. We fix 
the parameters of the effective quantum impurity Hamil- 
tonian using SU(2) and Galilean invariance, and find the 
exponents of threshold singularities of various dynamic 
responses [density and spin structure factors as well as 
A{k,u)]. 

Away from k = kp, the spectrum of the spinon mode 
e s (fc) departs from the linear one, v s (k—kp), and becomes 
a periodic function of k with period 2kp . The threshold 
for A(k,ui) is located at \u)k\ — e s (k~2nkp), with integer 
n such that \k — 2nkp\ < kp, see Fig. [TJ The spin-charge 
separation is preserved near the threshold at arbitrary k 
in the following sense: if the energy uj of, say, the ex- 
tracted fermion approaches the threshold Wfe, then the 
momentum of a created spinon is approaching k — kp. 
The rest of the energy, ~ \uj — uJk\ — > 0, is given to a 
holon; it inevitably resides near a Fermi point and may 
be described as a conventional linear LL. We express the 
exponents of various threshold singularities in terms of 
the derivatives of e s (k) with respect to k and p, the den- 
sity of the liquid [see Eq. (|11[) and Table Uj. The obtained 
exponents are valid at arbitrary k, including the Fermi 
points. 
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Near the Fermi points, the exponents for A(k, ui) ap- 
proach the universal values which depend only on K c . In 
the main region (n = 0), the exponent for A(k,oj < 0) 
coincides with the predictions of the linear LL (unlike in 
the spinless case [9(), while density and spin structure 
factor exponents approach 1/2. 

Let us first introduce the refcrmionization of the linear 
LL Hamiltonian. Its charge (c) and spin (s) parts sepa- 
rate from each other, Hq = H c + H s , and in terms of bo- 
son variables, these parts have the conventional form [l|, 



are given by 



2tt 



dx 



(1) 



for v = c,s. The canonically conjugate fields d x v 
and 4> v are the momentum and displacement operators 
of bosonic charge and spin density waves and satisfy 
\<j>v{x),d y ll {y)} = in5 l _ lu 5(x — y). The effects of the in- 
teraction are contained in the constants K v , but most 
importantly interaction yields a difference between the 
velocities of charge and spin modes, thus removing the 
degeneracy characteristic of the free-fermion system. For 
repulsive interactions at the SU (2) symmetric point, we 
can assume rather generally K s = 1 [H 0- Each of 
the Hamiltonians (c, s) may be refermionized in terms 
of free fermionic quasiparticles, the same way as it was 
suggested earlier [LJ, [lj] for spinless particles. One ob- 
tains a free linear Hamiltonian in terms of left- and right- 
moving fermionic spin and charge quasiparticles iba»(x) 
(a = R,L = ±), spinons and holons, respectively [14| : 



H = -i ^2 v u ^2 a J dx : ^ av {x)Vi) au {x) 

v=c,s a=R,L 



(2) 



The operators t^ av {x) have fermionic commutation rela- 
tions, and their vacuum state has unity occupation for 
negative (positive) momenta for a = R{L). The left- 
and right-moving components of the original fermions, 
defined by ipa-(x) = e~ lkFX ipLa + e. lkl ' x ipRa, can be ex- 
pressed in terms of refermionized quasiparticle operators 
as (a =t,l= ±) 



i> aa {x) oc 4) ac (x)F ac (x)ip a ° \x)F° s {x) 



(3) 



Here notations are such that ipasix) i s the creation (an- 
nihilation) operator of a spinon, and the string operator 
F as (x) has to be raised to power a. We did not write out 
the Klein factors explicitly, since they do not affect expec- 
tation values which define response functions. For exam- 
ple, the correlation function ({ip(x, t), ip '(0, 0)}) which 
defines A(k,u) conserves the spinon and holon numbers. 
The string operators 



F av {x) = exp <^ -ia / dy [6+ (y) + s 



are functions of the quasiparticle densities 
Puv{x) ='■ Vw^Vw^) - The effects of the in- 
teraction are contained in the phase shifts S± u , which 



5±u 



1 I 

2:r V2 ± 2 



8K V 



K s = 1 leads to <5_ s = and <5 +s /(2tt) = 1 - l/y/2. 
Within the linear spectrum approximation of the LL the- 
ory, the dynamics of the string operators is linear, and 
Eqs. (J2])-(j4]) lead to the conventional LL results. 

We first discuss the deviations from the LL results for 
the exponents of the spectral function at the charge mode 
in the vicinity of the Fermi point +kp. Similar to the 
spinless case [9(, one can describe the exponents by tak- 
ing into account only the quadratic nonlinearity in the 
spectrum of charge quasiparticles. Interactions between 
charge quasiparticles lead only to small corrections to 
the exponent at the charge mode for \k — Uf\ "C kp, 
and we obtain a singularity at the holon mass shell 
A(k,oj) cx [u> — e c (fc)]~ AIc , where 



2tt 



2tt 



5+c 
2tt 



(5) 



The nonlinearity in e c {k) changes the exponents com- 
pared to LL theory in an energy window of width ~ 
(k — kp) 2 around e c (k). Even at K c — > 1, the expo- 
nent /i c — > y/2-1, different fr om the LL prediction 1/2. 
Beyond this energy window, the usual LL behavior is 
recovered. 

On the other hand, interactions between spinons on 
the same branch are important in the description of 
the response functions at the spinon spectrum even for 
| k — kp\ <§; kp . The difference stems from the fact that 
the leading-order band curvature of the spinons is cubic 
[H| as the SU{2) symmetry enforces particle-hole sym- 
metry for the spinon mode. If one attempts now to treat 
the leading interactions between spinons perturbatively 
in the spirit of Ref . 0] , one finds that these interactions 
lead to ~ 0(1) changes of the exponents in the response 
functions, since the difference in the velocities of spinons 
is proportional to oc (k — kp) 2 unlike oc (k — kp) for 
holons. 

Nevertheless, it is still possible to describe the spectral 
function at the spinon mode using mobile impurity mod- 
els, as has been established in a number of articles 
We shall first apply the approach to calculate A(k,uj) for 
\k — kp\ <§; kp,uj < Oin the vicinity of the spinon mode 
v « e s (k). According to Eq. (j3|), for energies close to 
e s (fc), the configurations of lowest energy will contain a 
spinon hole at momentum kd = k — kp < and a holon 
at kp. Therefore, we can project the Hamiltonian onto a 
band structure which consists of a "deep" right-moving 
spinon hole d at momentum k^ as well as spinon and 
holon states, ip s and ip c , respectively, at the Fermi points. 
This leads to the Hamiltonian Hq + Hd + Hi nU where Hq 
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is given by Eq. ([2]), while the other terms read 



Hi, 



dx (P (x) [e s (fc) — ihvd\7]d(x), 
dx > V av (k)p av (x)d(x)d^ (x) 



a // 



(6) 



The sub-bands were linearized around kd and the Fermi 
points, and we used i/jr s oc ipR S + e lkd ' x d and ip av oc i/j al/ 
in all other cases. The velocity of the hole d is given by 
Vd — de s (k)/dk and V au (k) reflect the interaction of the 
hole with the modes near the Fermi points [16j . 

The interaction term (fiat can be removed using a 
unitary transformation [8l4l2j|. This will lead to addi- 
tional phase shifts A8 av = V av /{Vd — om v ) which de- 
termine the edge exponents. Except for A5 +s , these 
additional phase shifts are small since the correspond- 
ing V av vanish for k — > kp , while the denominator re- 
mains finite. The remaining phase shift A(5 +s can be 
fixed using the SU{2) symmetry [ltj. In particular, this 
symmetry requires identical exponents of the spin cor- 
relation functions S h (fc,w) and S zz (k,u), where e.g. 
S zz (k,co) = J dtdxe iUt - tkx (S z (x,t)S z (0,0)}. We can cal- 
culate the exponents of these functions at the spinon 
mass shell for general phase shifts 5+ s = 6 +s + A5+ s . 
The leading exponents for S ^(k,uj) and S zz (k,ui) read 

1 - [1/V2 =F 5* + J(2ir)] 2 , respectively. Hence, the SU{2) 
symmetry leads to 5*^ s = and rules out the interaction 
of the impurity with the low-energy spinons, and the LL 
result /x^ L (see Fig. Q} for the exponent of the spectral 
function remains valid for \k — kp\ <C kp. 

Even beyond the universal regime \k — kp\ <C kp , the 
exponents of the spectral function near its edge of sup- 
port survive, and may be determined using mobile impu- 
rity Hamiltonians. Above analysis in the limit \k — kp\ *C 
kp indicates that in the vicinity of the Fermi points the 
impurity has the quantum numbers of a spinon, and the 
impurity Hamiltonian is given by Eq. ([6]). The contin- 
uous evolution of edge exponents implies that the same 
Hamiltonian should describe the edge exponents even sig- 
nificantly away from Fermi points. As a consequence of 
such a spin-charge separation beyond low energies, for 
Galilean-invariant systems we can express all exponents 
in terms of the edge position only, as in the spinless case 
[lp| . For this purpose, we rewrite Eq. as 



H tnt = I dx 



V R V- 



2tt 



V L V- 



2tt 



d(x)d*(x), 



since we saw previously that the coupling constant to the 
spin sector vanishes. The argument was based on SU(2)- 
symmetry and is valid beyond the low-energy regime. 

The interaction term is removed [lfj by a uni- 
tary transformation U'(H C + Hd + Hi nt )U , where 



U 



exp{i J dx[ 



2ir 
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FIG. 1: (Color online) Structure of the spectral function 
A(k,oj) in the (fc,w)-plane (A) and along a cut for fixed 
< k < kF (B,C). (A) Shaded areas indicate the regions 
where A(k, u>) is nonzero. (B) Luttinger liquid (LL) results: in 
the vicinity of the spinon mass shell, A(k,ui) oc {u — v B k)~^- 
where p^ = § - \ (K c + K' 1 - 2) . (C) Schematic behavior 
of the spectral function beyond the LL approximation away 
from the Fermi points: the exponent po,- is different from 
P- L , the charge mode is smeared out, and for uj > the re- 
gion of finite support starts from — e s (k) with an exponent 
/io,+ instead of from — e c (k). 



6> C V K c )]d(x)d^ (x)} with phases A8± c defined by 

(V L T V R )K^ 2 = -A6- C (v d + v c ) ± A5 +C (v d - v c ). 

The fermionic operator can be written as tp^ = 
^ e i(e c -ip c )/^ ( w j|; n t ne second factor coming from a 
holon at k — ¥ kp), and then its correlations can be eval- 
uated using the operator U similar to the spinless case. 
The phases A<5-|- c fix the exponents of the spectral func- 
tion as well as of the dynamic density and spin structure 
factors [1411. 

In order to relate Vl,r to the spectrum e s (k), we first 
consider the shift in energy due to a change in density. 
A uniform variation of the density by 5p corresponds to 
a finite expectation value (V(j) c ) = —ir5p/^/2, and leads 
to the shift in the single-particle energy 



Se s (k) 



de s (k) 
dp 



dp 

dp" 



Sp, 



(7) 



ir(: 



/7v 



where p, denotes the chemical potential. We now need to 
calculate the same change in energy using Eq. ©. Un- 
like the spinless case [loj , the spinon momentum changes 
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A(k,u ^ 0) 


(2n - l)k F < k < (2n + l)k F 


P>n,± 




' (2n + iyj£ ! ^ 
x/2 


+ s*\ 

2tv J 
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f 1 5+ - S-\ 
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S(k,u) 


2nk F < k < 2(n + l)k F 


u DSF 


1 i 

2 ~ 2 


\ V2 2, J 


*~\ 


(8 s ;- 

V 2tt 
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TABLE I: Exponents for the spectral function A(k, uj) (see Fig. Q] for notations) and the dynamic density structure factor 
S(k,u) at the edge of support. The exponents are determined in terms of the phase shifts 5± = AS± c (k — 2nk F ) and 
5± = A8± c [(2n + l)k F — k] calculated in Eq. (j 1 1 p . Moreover, m± = (n + l/2± 1/2) mod 2. The exponents for the spin structure 
factor S zz (k,uu) coincide with the exponents for S(k,w). 



under variation of the density, since the total momentum 
is fixed while the holon momentum changes following the 
shift of the Fermi point. Calculating the shift in the en- 
ergy of the spinon d as well as the shift in energy of the 
holon at the Fermi point, and comparing with Eq. ([7]). 
one finds 



V R + V L _ de s (k) , 7Tde s (k) 



2-y/2 



dp 



+ 



2 dk 



(8) 



The second relation can be derived using Galilean in- 
variance, which predicts flp| that a uniform change in 
velocity u should lead to a change in energy of 



Se s (k) 



mu 



k 



de s (k) 



Ok 



(9) 



where m is the bare mass. On the other hand, this 
change in velocity u will lead to a finite expectation 
value, (V0 C ) = V2mu. Because of the shift of the Fermi 
point, this leads to an energy shift of a holon at the 
right Fermi point K c v c mu = vftuu, and a spinon en- 
ergy shift —muvd due to the change of the spinon mo- 
mentum. In addition, the interaction Hamiltonian yields 
a shift y/2mu{VL — Vr)/(27t). Combining these terms 
leads to 



Vl-Vr = k_ 
V2tt 



(10) 



Equations ([8]) and ((10]) allow us to express Vl,r in terms 
of the derivatives of the single-particle spectrum e s (fc). 
For the phase shifts, we thus obtain the result 



AJ± c (fc) 
2tt 



± 



k — k £r 
m-JW, 



(» 



de s (k) 
dp 



de s {k) 
dk 



2^2 (2%M T J^- 

v V 7r ok 1 mK c 



(ii) 



The predictions (| 1 1 1) and the absence of the coupling of 
the spinon impurity to the low-energy spinons can be 
explicitly checked for the case of the integrable Yang- 
Gaudin model based on its finite size spectrum [17j . 
Moreover, for k -> k F , Eqs. © and (JTOJ) yield V L = 
Vr — which reflects the absence of interactions between 
holons and spinons within the linear LL theory. 

Finally, let us present the exponents of the spectral 
function in the regions (2n — l)kp < k < (2n + l)kp for 
integer n. For uj < 0, in the vicinity of e s {k), the sym- 
metry of the edge position in k leads to a spinon hole at 



momentum kd, n — k — (2n + 1)/cf < 0, a holon at kp as 
well as additional excitations which absorb the remaining 
momentum 2nkp. In terms of the original fermions, these 
"umklapp" excitations contain particles or holes at the 
Fermi points, and are thus characterized by four param- 
eters. Fixing the total momentum at 2nkp and requir- 
ing zero total charge and spin leaves one parameter m_, 
where the umklapp contribution to ipa- can be represented 
as {^ Rf ^L-ty {n+m - )/2 {^{ l ^Ri) {n ' m - )/2 , and thus m_ 
has to satisfy the selection rule m_ = n (mod 2). The 
exponent of the spectral function is now determined by 
the phase shifts (fTTj) taken at momentum k — 2nfcF.[10l| 
For ui > 0, one can also describe the exponents using the 
impurity Hamiltonian with the same parameters, but the 
states which determine the exponents are different. Fol- 
lowing the same line of arguments as previously, one finds 
exponents which are formally identical to the uj < case, 
but where m + now has to satisfy a different selection rule, 
TYi-if. = (n+ 1) (mod 2). The leading exponents stem from 
m± with smallest absolute value allowed by the selection 
rules, and the results are shown in Table HI 

The exponents of the spin structure fac- 
tor S zz (k,uj) and the dynamic structure factor 
S(k,u) = Jdtdx e i( - wt - kx ) (p(x,t)p(0,0)) can also 
be calculated similarly. It turns out that the exponents 
for S(k,uj) coincide with those for S zz (k,uj), and they 
are shown in Table Q] 

In conclusion, we have calculated the exponents of dy- 
namical correlation functions of ID spinful interacting 
fermionic systems beyond the approximation of a linear 
spectrum, shedding light on the fate of the spin-charge 
separation away from Fermi points. In the low-energy 
sector near the Fermi points, we found universal phase 
shifts which control the exponents and depend only on 
the LL parameter K c . Beyond the low-energy regime, 
we were able to establish phenomenological relations be- 
tween phase shifts and properties of the spinon spectrum. 

We thank F. Essler and R.G. Pereira for useful discus- 
sions. We acknowledge support by the NSF DMR Grant 
No. 0906498, the Nanosciences Foundation at Grenoble, 
France, and the Swiss NSF. 

Note Added. - Recently, the preprint 18] appeared, 
where the exponents for the density and spin structure 
factors for k <C kp were derived. 
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